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Abstract. We study the quantum dynamics resulting from preparing a one- 
dimensional quantum system in the ground state of initially two decoupled parts 
which are then joined together (local quench). Specifically we focus on the transverse 
Ising chain and compute the time-dependence of the magnetization profile, mi{t), and 
correlation functions at the critical point, in the ferromagnetically ordered phase and 
in the paramagnetic phase. At the critical point we find finite size scaling forms for 
the noncquilibrium magnetization and compare predictions of conformal field theory 
with our numerical results. In the ferromagnetic phase the magnetization profiles are 
well matched by our predictions from a quasi-classical calculation. 
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1. Introduction 

Recently there has been an increased interest in nonequilibrium relaxation processes in 
closed quantum systems following a sudden change of the parameters of the Hamiltonian 
(quantum quench) P[2l[3lli[5l[6l[7l[I7lElElEQlEl[ia[l3lE3Ea[l6l[ig 
Theoretically during the initial period (t < 0) the system is described by a Hamiltonian 
"Ho with ground state I^I/q), which is suddenly changed to a new Hamiltonian, "H, for 
time t > 0. The new state of the system - in the Schrodinger picture - is time-dependent: 

|vl>(t))=exp(-^t-H)|vl>o), (1) 

and so is with an operator, say O, which is expressed at time t > in the Heisenberg 
picture as: 

0{t) = exp{-itn)0 exp{itn) . (2) 
Generally one is interested in the relaxation of the local order-parameter, 

m,(t) = (^o|a(t)|^o), (3) 
at site, r or the time- dependence of the autocorrelation function: 

a(t,t') = (^oia(t)a(Oi^o). (4) 

After long times, t ^ 1, the system is expected to relax to a stationary state, in which 
one measures the correlation function: 

a(ri,r2) = (vl/o|a,(t)a,(t)|^o). (5) 

Most often one considers global quenches, when the parameters are modified uniformly 
at all points of the sample. Experimentally this process can be realized with ultracold 
atomic gases in optical lattices |2T | |22 | [23 | I24]. After a global quench a quantum system is 
expected to relax to a thermal (or thermal-like) state, such that the local magnetization 
(as well as the autocorrelation function) vanishes exponentially, mr{t) ~ exp(— t/r), 
with a relaxation (or phase coherence) time r. Similarly, the correlation function in the 
stationary state behaves as: 

Ct{ri,r2) ~ exp(-|ri - r2|/0 , (6) 

with a finite non-equilibrium correlation length, ^. This thermalization of the system 
is generally explained in terms of quasi-particles [H [25], which are emitted during 
the quench at each point of the sample and which travel with a constant velocity, v. 
Quasi-particles which are originated at nearby sites (which are within the correlation 
length) are quantum entangled, but other quasi-particles are incoherent. Incoherent 
particles reaching a reference point, r, will result in the reduction of the value of 
a local observable and this process is responsible for the exponential decay of the 
(auto) correlation function. In a system with boundaries these quasi-particles are 
reflected at the boundaries, which results in more complicated time-dependence, among 
others reconstruction of the local magnetization [19j . 
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Besides global quenches one also considers local quenches, when the parameters are 
modified only at a restricted region of the system [26| [27] . Experimentally this type 
of problem can be realized in the x-ray absorption problem in metals [29j, where the 
creation of a hole plays the role of a local defect and when a conduction electron fills 
the hole this potential is suddenly switched off. Nonequilibrium dynamics following a 
local quench has been first studied in the context of the entanglement entropy, S, in 
critical quantum spin chains [26| [27 t |28 ] . If we consider two disconnected half chains for 
time t < 0, which are suddenly joint together for t > 0, then the time-evolution of the 
entanglement entropy of the two half chains is found to evolve in time asymptotically 
as: 



where c is the central charge of the Virasoro algebra and Ci is a non-universal constant. 
In view of Eq.Q basic informations about the critical properties of quantum spin chains, 
such as the value of the central charge, can be obtained through nonequilibrium local 
quench dynamics. The result in Eq.(I7]), which has been first obtained numerically for 
a free-fermion model |26], has been derived through conformal invariance |27]. It has 
also been generalized for different positions of the defect [261 [27] . as well as for varying 
strength of the defect in quantum Ising and XX-chains [30]. For a large finite chain of 
total length, L, the appropriate expression is conjectured to be [30] : 



The conformal mapping, which has been used for the calculation of the entanglement 
entropy can be applied [27] to study the time- dependence of one-point functions (such as 
the local magnetization in Eq.Q), as well as two-point functions (such as the correlation 
function in Eq.(j5])) in a critical quantum spin chain. In contrary to global quenches 
these results indicate a power-low relaxation as well as power-low type correlations in 
the stationary state. This is understandable within the quasi-particle picture, since at 
a local quench the quasi-particles are expected to be emitted only at a restricted region 
of the sample and therefore these are quantum entangled. If these quasi-particles arrive 
after time t to different positions of the chain, say ri and r2, will result in correlations. 

Besides the conformal results which are mentioned above several problems 
connected with local quench dynamics are still unexplored. To our knowledge there 
are no numerical investigations, which could be used to check the validity of the 
conformal conjectures. There are no results about the autocorrelation function and 
no information is known about nonequilibrium relaxation following a local quench to 
the ordered or to the disordered phase of a quantum spin chain. In this paper we 
aim to fill this gap and perform detailed numerical investigations about these open 
questions. As a model we use the transverse Ising spin chain, which is a prototypical 
quantum system having ordered and disordered phases as well as a quantum critical 
point [31j. Using free-fermionic techniques [32t [33| |3T] we have studied numerically the 
nonequilibrium relaxation processes for the magnetization, the autocorrelation and the 



5(t) = -logt + ci, 



(7) 



^L(t) = - log — sin — + ci . 



(8) 



Non- equilibrium quantum dynamics after local quenches 



4 



equal-time correlation function in large finite systems. Most of our investigations are 
performed at the quantum critical point, so that we could compare our results with the 
conformal predictions. However, we also studied quenches to the ordered and to the 
disordered phases and these results are then compared with semi-classical calculations.. 

The structure of the paper is the following. The model and the numerical method 
of the calculation is introduced in Sec|2l Results of the calculations at the critical point, 
in the ferromagnetic and in the paramagnetic phases are presented in Sections [3l H] and 
El respectively. Our conclusions are summarized in Secj6l 



2. Model 



The system we consider in this paper is the transverse Ising chain of finite length L with 
open boundaries defined by the Hamiltonian: 



1=1 1=1 



(9) 



in terms of the Pauli- matrices ' at site / and bond strengths J;, which are all equal 
Ji = 1 except the central bond Jl/2, which is 7^/2 = for time t < 0. Similarly, 
the transverse fields are homogeneous, hi = h, except at the central sites, which are 
hL/2 = = for time t < for fixed-spin boundary conditions. (For free-spin 

boundary conditions these are hL/2 = ^l/2+i = h for time t < 0, too.) Generally we 
measure distances from the defect and use the variable: 

r = /-§. (10) 

The order-parameter operator of the system is, erf, what should be inserted in the 
formulae of the autocorrelation and the correlation functions, see in Eqs.(jl]) and 
([5]), respectively. For the local magnetization in Eq.([3]) one should have: mr{t) = 
lim6_j,o+ b{^oWr{t)\'^o)b, where 1^^0)6 is the ground state of the initial Hamiltonian ([9]) 
in the presence of an external longitudinal field b. According to [34j this can be written 
as the off-diagonal matrix-element of the Hamiltonian ([9]): 

m,(t) = (vl/o|<(t)|vl/i), (11) 

where |\E'i) is the first excited state of the initial Hamiltonian (t < 0). For fixed-spin 
boundary condition the ground-state is exactly degenerate with |\E'i). In the initial state 
and in the thermodynamic limit, L — )■ oo, there is spontaneous ferromagnetic order in 
the system, mj.(0) > 0, for h < he = 1. On the contrary, for stronger transverse fields, 
h > he, the magnetization is vanishing. At the quantum critical point the magnetization 
vanishes as a power-law: 771^(0) ~ L~^, for bulk spins: r = 0(1), and r«±L/2(0) ~ L~^% 
for surface spins. Here the magnetization exponent x and the surface magnetization 
exponent Xg is known exactly: x = 1/8 and Xg = 1/2. 

The Hamiltonian in Eq. ([9]) can be expressed in terms of free fermions [32| |33| |3T] , 
which has been used in previous studies of its non-equilibrium properties [35| [T8l [T9] . 
The magnetization profile, as well as the (auto) correlation functions can be expressed 
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in terms of Pfaffians, which are then evaluated as the square-root of an antisymmetric 
matrix, which has a rank 0{L). In the following Sections we use these techniques to 
calculate different nonequilibrium quantities. 

3. Local quenches at the critical point 

3.1. Conformal field theory - a reminder 

Here we recapitulate the basic results by Calabrese and Cardy |27j about the use of 
conformal field theory for local quenches at the critical point. The system is represented 
in a space-time region with coordinates: (r, t) and the two parts of the chain are 
decoupled for t < and these are joined at t = and we measure the properties 
of the system at t > 0. In the path integral formalism one introduces damping 
factors: exp(— e'H), so that in Eq.(l2]) we have expl—tfH) exp(— zfH — eH) and 
exp(«t'H) — )• expltfH — e'H). For computational simplicity the operators are inserted 
at imaginary times: r = it and one works in the complex plane z = r + it. Here, 
due to the local quench we have two cuts starting at ±ze and ending at ±2oo. This is 
represented in the left panel of Fig|T] 



2e 



w 



Figure 1. (Color online) Space-time region z = r + it with two cuts starting at ±«e 
(left) which is mapped to the half-plane w through the conformal transformation in 
Eq. (fT2|) . At the end of the calculation one should take r it. 



This geometry is mapped by the Joukowsky transformation: 

into the half plane. Re w > 0. At the boundaries of the slits in the z-plane as well as at 
the surface of the w-plane one should impose boundary conditions, which can be either 
free or fixed spin boundary conditions. In the w-plane the asymptotic form of the one- 
and two-point functions are generally known due to conformal invariance [36 | [ST f [38 | HQ] . 
These results are then transformed back to the z plane and at the end of the calculation 
one analytically continues the final result as r — it. 
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We note that the conformal results obtained in this way are vahd for infinitely long 
chains, i.e. for L — t- oo. 

3.2. Magnetization 

We have calculated the relaxation of the magnetization in the transverse-field Ising 
chain following a local quench starting with two different type of initial state. First, we 
consider an initial state with fixed spins at the defect, in which case the initial defect 
magnetization stays finite. In this type of setting one can make a direct comparison 
with the conjectures of conformal invariance. In the second type of calculation in the 
initial state we use free boundary conditions at the defect. Then we calculate the time- 
dependence of the off-diagonal order-parameter in Eq.f llip . which has a vanishing initial 
value at the defect in the thermodynamical limit. 

3.2.1. Fixed-spin boundary condition In this section in the initial state, t < 0, the two 
half chains are prepared with fixed spins at the boundaries: (cr^/g) ~ ('^l/2+i) ~ 
which is obtained by fixing /il/2 = ^l/2+i = in the Hamiltonian in Eq.(l9]). We use the 
superscript, to refer for fixed-spin boundary condition. The magnetization profile 
in the initial state is known from conformal invariance [37t [38] : 



which behaves close to the defect: m[^\t = 0) = y4|r|~^', |r| <^ L/2. This is the 
well-known result by Fisher and de Gennes |39) . 

After the local quench, for t > 0, we take = ^l/2+i = 1 (and Jl/2 = 1) and 
study the evolution of the magnetization. In the limit |r| ^ L/2 the conformal mapping 
in the Sec lS.ll leads to the result |27] : 



This can be interpreted, that for t < r the magnetization keeps its initial value until the 
quasiparticles from the defect arrive at t = r and afterwards for r <^t <^ L/2 the decay 
is given by m^^\t) ~ t''^^ . This decay involves twice of the magnetization exponent 
and this behaviour is similar to that of the equilibrium autocorrelation function at the 
critical point. 

In order to check the conformal result in Eq. (fT^ we have calculated the time- 
dependence of the magnetization in finite chains of length up to L = 256. For the 
largest chain and for different values of r the relaxation of the magnetization is shown 
in Figj2l In agreement with the quasi-particle picture and with the results of conformal 
invariance the local magnetization stays unchanged until t < r, which is followed by 
a fast decrease. Due to the finite size of the system the magnetization has a finite, 
L-dependent limiting value and for t > L/2 the magnetization starts to increase. 

We have studied in more detail the behavior of the magnetization at the defect: 
/ = L/2, i.e. for r = 1. For short times, 1 ^ t ^ L/2, the decay is compatible 




(13) 




t <r, 
t > r . 



(14) 
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Figure 2. (Color online) Temporal evolution of the magnetization at the critical 
point after a local quench having fixed spins at the defect in the initial state in a finite 
system of length L = 256. 




Figure 3. (Color online) Temporal evolution of the magnetization at the defect after 
a local quench from a fixed spin at the critical point. In the main panel the scaled 
profiles are shown for different values of L, where the curve with the dashed line is 
given in Eg. p^ . In the inset the log- log plot of m^^^t, L) is shown for L — 256, where 
the conformal prediction for L — >■ cx) is given by the dashed straight line. 



with the conformal prediction: rui (t) ~ t^^/^ is seen in the inset of Figj3] in 
which the magnetization is shown as a function of time in a log- log plot. For 
longer time the boundaries of the chain start to influence the relaxation and the 
critical defect magnetization is expected to satisfy the scaling behavior: m^^\t,L) = 
b''^^mi^\t/b^ , L/b), when lengths are rescaled by a factor 6 > 1. In this relation the 
dynamical exponent is 2; = 1 and by taking the scaling factor as 6 = t in the limit L S> t 
we recover the conformal result: m^^-*(t) ~ t~'^^. It is more interesting to take b = L, 
when one obtains: mi^\t,L) = L'"^^ ^^^\t / L) . Here the scaling function /ii^''(r) for 
small T = t/L behaves as: /iQ^'*(r) ~ r^^^. Numerical results for the scaling function 
for different values of L are shown in the main panel of Fig|3l which is found to be well 
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Figure 4. (Color online) Temporal evolution of the local magnetization profiles at 
the critical point after a local quench with free boundary conditions in the initial state. 



approximated by the function B [sin(7rr)] thus we have the conjecture: 



-2x 



sm TT 



L 



-2x 



(15) 



3.2.2. Free-spin boundary condition In this section in the initial state the two half 
chains have free boundary conditions, which means that in Eq.(l9]) we have hL/2 = 
= 1 and Jl/2 = 0. In this settings no conformal results are available, therefore 
we study numerically the properties of the profiles of the off-diagonal order parameter, 
mr{t) in Eq. lfTT]) in finite systems. The results are shown for L = 128 in Fig. HI 

Initially (at t = 0) the ground state magnetization profiles are identical and 
independent in both disconnected parts of the system. The functional form is known 
from conformal invariance [40j and the complete profile has the finite size scaling form 

sm . (loj 



mr{t = 0, L) oc L~ 



Close to the defect, i.e. for \r\ <^ L/2 the behavior of the profile follows from the scaling 
relation, that mr{t = 0,L) = b~^mr/b{t = 0,L/b), with a rescaling factor, b > 1. Now 
taking b = L, we arrive to the form: mr{t = 0,L) = L~^m{r / L) , where the scaling 
function, m(p) for small argument behaves as m(p) ~ p^""^ . This is compatible with 
the conformal result in Eq. (fT6l) . if we replace p by its sinusoidal extension: sin{27[p). 

In Fig. m one observes a quasi-periodic time- dependence of the profile with the 
characteristic initial double peak being exchanged against a single peak at times 
T = L/2, 3L/2, . . . and re-occurring at times T = L, 2L, . . .. Let us focus first on the 
time- dependence of the magnetization at the central cite, mi{t, L), which is expected to 
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Figure 5. (Color online) Scaling plot of the local magnetization at the defect, ma{t), 
after a local quench at criticality for different system size L. The conjectured result in 
Eq. (fT7|) is shown by a full line. Inset: Log-log plot of mo{t) vs. t for L = 256. The 
straight line is proportional to t^^^. 



satisfy the same type of scaling relation as m^\t, L), thus: mi{t, L) = b''^^mi{t/b, L/b). 
As before taking b = L we arrive to: mi{t,L) = L~'^^]li{t/L), where the scaling 
function, /ii(r), for small argument should behave as fliir) ~ ^xs-2x^ ^j^jg 
we obtain: mi{t,L) = L~^=t^=~^^', which is in agreement with the L-dependence of the 
magnetization at the defect at t = 0, see in Eq.(fT6l). Furthermore we have for the 
relaxation for small times: mi(t) ~ t^/'*, which agrees well with the numerical data 
shown in the inset of Fig. [5l The form of the scaling function can be conjectured using 
the same substitution, p — )■ sin{2TTp), as for t = 0. In this way we obtain: 



mi{t, L) oc L 



-2x 



27rt ^=-2^ 



sm — — 
Ij 



(17) 



Fig. ([5]) displays a corresponding finite size scaling plot for the first period which shows 
a good data collapse (a corresponding scaling plot for larger values of t/L is equally 
good, data not shown) and thus confirms our conjecture (fT7|) . which probably can be 
derived rigorously from conformal invariance. 

Next we consider t = L/2, when the profile has its maximum at r = and it is 
minimal at the two ends of the chain, r = ±L/2, i.e. at Z = 1 and I = L. Let us 
consider the profile for small / and use the scaling transformation: mi{t/L = 1/2, L) = 
b-^^'mi/bit/L = 1/2, L/b), which leads to the result: mi{t/L = 1/2, L) = L-^'^pil/L) , 
with b = L. The scaling function, p{y), for small argument is expected to behave in 
the same way as for t = 0, thus fi{y) ~ y^^-x _ Furthermore having the substitution: 
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y — sin(27r|/) we arrive to the conjecture: 
rurit = L/2,L) oc L'^^ 



vrr 
cos — 
ij 



where we use the variable r. Fig. [6] displays a corresponding finite size scaling plot 
which shows a good data collapse and thus supports our conjecture fll8l) . 

Combining the scaling forms ( IT7|) and ( !T8|) the ratio i? = '^/=j;^/2('^ ~ 't)/^\=t{^ = 
is given by 

= fsinvr^) . (19) 



L 

3.3. Spatial correlations 

The spatial correlation function, Ct{ri, in Eq.([5]), has also been studied by conformal 
field theory |27| and various analytical predictions have been made in an infinite system 
in the continuum limit. In this section we will compare our results for finite lattice 
systems with free boundaries with these predictions. We note that before the quench 
the two halves of the system are disconnected and free, i.e. hL/2 = /iL/2+1 = 1 and 
Jl/2 = 0. Without the restriction of generality we take ri > \r2\ > and consider two 
cases: 1) both sites of reference are at the same side of the defect (r2 > 0), and 2) the 
two sites are at different side of the defect (r2 < 0). 

Case 1: r2 > For short times, t < r2(< ri), the behavior of the correlation 
function can be obtained in the quasi-particle picture. In this case the quasi-particles 
(which can be called as signals) starting at the defect and propagating with a speed 
V = 1 does not reach any point of reference, thus the correlations keep their initial 
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value: Ct{ri,r2) = Co{ri,r2) oc |ri — r2|~^^'. This result follows also from conformal field 
theory and consistent with our numerical data depicted in Fig. [7] a and b. 

For intermediate times, r2 < t < ri, the signals reach the closest site at r2 but not 
yet the remote site at ri, see [7] a and b. The prediction from conformal field theory [2] 
is 

f (ri + r2)(r2 + t)e 



Ct{ri,r2) 



(ri - r2)(ri - t)4ri(t2 
2ri(r2 + t) 
(ri + r2)(ri + t)^ 




X F (^Jri^±I>] (20) 



F(l)={Jl + Vn-Jl-V^l)/-^~Vn forij^O (21) 



with 



and e is the regularization parameter in Eq. fll2p . In Fig. [7b a comparison with our 
numerical data for a specific vi and r2 is shown. As expected for a lattice model our 
data display characteristic oscillations in the considered regime around the monotonous 
continuum result (120|) . 

For a more quantitative comparison we first look at the asymptotic behavior of 
Ct{ri, r2) for fixed small t and r2- The conformal field theory predicts according to (120|) 
and (EH) 

r2) oc rf 2^-^^ = r^"^/^ for t < r2 < n (22) 

which agrees well with our lattice result as shown in Fig. [T^. Moreover, for fixing the 
ratios P2 = r2/t < 1 and pi = ri/t > 1 ( |20ll2T]) predicts in the limit t — )■ cxd 

a(ri, r2) oc t"^^' = t"^/^ for p2 < 1 Pi > 1 t ^ oo . (23) 

This is also consistent with our finite lattice data as seen in Fig. [T]! for p2 = 1/2 and 
Pi = 2. 

In the long-time regime, for t > ri the signals have also reached the remote site 
at ri, and both points of reference has the information of the joining of the two halves 
of the system. Consequently Ct{ri,r2) approaches a time-independent value (c.f. Fig. 
I20k and b), given by the equilibrium bulk correlation function, which is proportional to 
\fi — r2\ . 

Case 2: r2 < 0, i.e. the two reference sites are at the different sides of the defect. 
Initially, at t = 0, the two magnetization operators are located in separate, independent 
subsystems, thus their correlations should vanish. As long as after the quench no signal 
reaches the site at r2 (i.e. for t < \r2\) the correlation function should be expected to 
be zero. (Strictly speaking Ct{ri,r2) can be exponentially small, the only contributions 
coming from signals propagating outside the "light-cone".) This is clearly visible in our 
results shown in Fig. [Hti. Surprisingly this disagrees with the conformal field theory 
result [27] which predicts 

G(ri,r2) = -^-Ff^^4Sl^^^^l fort<|r2|. (24) 



4rir2|^ \{rl-t^){r^-t^) 



Non- equilibrium quantum dynamics after local quenches 



12 



B 



0.3 



0.2 



0.1 



-1 
-1.5 

-2 
-2.5 

-3 
-3.5 

-4 



1.5 



20 







4, r.|=24 

r2=10, ri=30 

r2=20, ri=40 


1 \ ^J\: \i 



40 



60 



80 



O 



100 



D 





-a/4 




t=4 




t-5 




t=6 


x.:;;Cx '-1=2, 


t=7 



2.5 



3 3.5 
log 



4 4.5 



5.5 




10 20 30 40 50 60 70 80 
t 




Figure 7. (Color online) A: Spatial correlation function Ct{ri,r2) as a function of 
time t for both sites of reference being at the same side of the defect (i.e. ri > r2 > 0). 
The system size is L — 256, as in the other panels. B: Comparison of the data of A 
for r2 = 10, ri = 30 with the conformal field theory result ([20)) . C: Ct{ri,r2) in the 
intermediate time regime (r2 < < < ri) for small r2 and i as a function of ri. In the 
limit ri — > oo this should approach according to the conformal field theory result 

([22| . D: C{ri,r2;t) in the intermediate time regime for r2 — t/2 and ri = 2i as a 
function of t that should decay as according to the conformal field theory result 

(El. 



This could be rectified by setting e = 0, but this would be at variance with the results for 
r2 > 0, for which the regularization parameter should be non-vanishing. The problem 
with the conformal derivation could be related to the fact, that in the transformation 
f fT2|) points in the z plane with Re(z2) = r2 < and \r2\ ^ e are transformed to \w2\ <^ e. 
However, the scaling form of the correlation function in the semi-infinite w-plane is valid 
only in the continuum limit, i.e. for R.e{w2) ^ 1. 

After the signals reach the site at r2, i.e. for |r2| < t < ri, the conformal field theory 
prediction is [27J 

/ (ri + r2)(r2 + t)e V 



Ct{ri,r2) 



(ri — r2)(ri — t)4ri(t^ — rf) 



^f(-^^^^^±^] for|r2|<t<ri. (25) 
\in + r2)iri + t) J 

In Fig. [S)3 a comparison with our numerical data for a specific Vi and r2 is shown. As 
expected for a lattice model our data display characteristic oscillations in the considered 
regime around the monotonous continuum result (1251) . 
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Figure 8. (Color online) A: Spatial correlation function Ct{ri,r2) as a function of 
time t for the two sites of reference being at different sides of the defect: ri > 0, r2 < 
(''I > l''2|)- The system size is L = 256, as in the other panels. B: Comparison of 
Ct{ri,r2) with the conformal field theory result ([25]) . C: Ct{ri,r2) for small \r2\ and t 
as a function of ri. In the limit ri — > cx) this should approach according to the 

conformal field theory result ((22|). D: Ct{ri,r2) for r2 = t/2 and ri = 2i as a function 
of t that should decay as t^^^^ according to the conformal field theory result (|23|). 



For fixed small t and r2 the conformal field theory predicts according to (l25l) and 
f l2T]) again for Ct{ri, the asymptotic ri-behavior (122|) as in the case r2 > 0. This agrees 
well with our lattice result as shown in Fig. [8b. Moreover, for fixed p2 = r2/t > —1 and 
Pi = Ti/t > 1 ( !25|) predicts in the limit t — )■ oo the same asymptotic t-behavior ( 123|) as 
in the case r2 > 0. This is also in agreement with our lattice result as shown in Fig. |8jl. 
Finally for t > vi the signal has also reached the site at ri and Ct{ri,r2) approaches 
again the time-independent equilibrium bulk value (c.f. Fig. [8]), which is proportional 
to |ri — r2|~^^. 

Summarizing our results confirm the predictions of the conformal field theory ^27] 
for the spatial correlation after a local quench at the critical point, except for the case 
r2 < 0, t < |r2|. 

3.4- Autocorrelations 

We have also calculated the autocorrelation function, Gr{t,0), as defined in Eq.(jll). In 
the following for simplicity we shall omit the second argument and use the notation 
Gr{t). Using the quasi-particle picture we have the following expectations. Before the 
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signal reaches the reference point, t < r, the autocorrelation function is the same as in 
the equilibrium bulk system, thus asymptotically Gr{t) oc t-^x/z _ ^-i/4^ p^j, t > 
when the signal has passed the reference point the equilibrium bulk decay of Gr{t) 
should continue, thus in the complete time-window this behavior is observable. 

0.15 



0.13 
0.11 



0.09 

50 100 150 200 

t 

Figure 9. (Color online) Autocorrelation function Gr{t) after a critical quench for 
different system sizes. The straight line corresponds to the equilibrium power law 

Our lattice results in Fig. [9] are in agreement with these expectations. In a finite 
lattice of length L the autocorrelation function has a minimal value of ~ L~^/^. 

4. Quench in the ferromagnetic phase 

In equilibrium in the ferromagnetic (FM) phase {h < 1) there is spontaneous order and 
the bulk rUb and the surface, m^, magnetizations are given by: 

m, = (l-/i2)-, = (1 - /i^)-^-, (26) 

respectively (the correlation length exponent is z/ = 1). The magnetization profile has an 
exponential variation in the surface region, the size of which is given by: C,s ~ 
close to the critical point. Here we follow the same protocol as in Sec J3.2.2l for t < we 
cut the system into two halves with free boundary conditions, which are then (at t = 0) 
joined together with the coupling = J = 1- We measure the temporal evolution 
of the magnetization profile, mr{t), as defined in Eq. (fTT]) . The finite lattice results are 
depicted in Fig. [inifor h = 0.5 and L = 128. In the initial state the magnetization profile 
{mr{t = 0)) is essentially constant and given by m^, except close to the boundaries and 
to the center. For t > one observes again a quasi-periodic pattern and the period of 
time T{h) is found to increase with increasing L and decreasing field h. It turns out 
that the spatio-temporal evolution of the profile can be understood even quantitatively 
within a quasi-classical picture. 

As elaborated in |41| within the FM phase the dynamics of the local order parameter 
in a system with boundaries after a global quench is very well described by a quasi- 
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Figure 10. (Color online) Temporal evolution of the local magnetization profiles 
after a local quench within the ordered (FM) phase, here h = 0.5, L = 128. 



particle (QP) picture, where QPs are kinks that move with velocity ±Vp through the 
system. The QP energy is 



'l + h"^ -2hcos{p) (27) 
with p = {2n — 1)tt/2L {n = 1, . . . ,L) and the QP velocity is 

.,= 1^ = ^^. (28) 

dp ep 



The maximum velocity is given by: 

.2 



V 



1 + /l2)2 + 12/l2 -{l + h' 



/6 , (29) 



for small h it is Vj^ax ~ h- 

Consider now a QP, or kink, with momentum p. It moves uniformly with velocity 
Vp until it reaches one of the boundaries, where it is reflected and moves with velocity 
—Vp thereafter, and so forth. The trajectory of the kink is periodic in time, after a time 
2Tp with 

Tp = L/vp (30) 

(including a reflection at the right and left boundary) it returns to the starting point 
Xo with the initial direction and velocity Vp. After a global quench QPs emerge pairwise 
at random positions Xq with velocities +Vp and —Vp [T8| E] and therefore we assume 
that after the local quench studied here QPs emerge also pairwise, but exclusively at 
the central site, where the defect is located before the quench. 
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The time dependent decay of the local magnetization mj.{t) at a position r is then 
determined by the probability q with which any given kink trajectory passes until time 
t the site r an odd (!) number of times: Each passage of one of the two trajectories flips 
the spin at site r and an even number of flips does not change the magnetization of site 
r. Once q{t) is known the magnetization is given by 

mr{t) = nirit = 0) ■ exp(-2g(t)) . (31) 

The probability q{t) is expressed as 

q{t) = - r dpU{h)q,{t) (32) 

TT JO 

where qp{t) is the probability that any given QP trajectory passes the site r an odd 
number of times, and fp{h) is the probability with which QPs with momentum p are 
generated (per site), and we take assume that it is identical to fp after a global quench. 
This is given for small h as [41j : 

fp{h) = -h's\n\p) (33) 

To calculate one concentrates first on times t < 2Tp and on lattice site r < 
- the whole profile is of course symmetric with respect to a reflection at the center 
mr{t) = m^r{t)- The QPs emerging at the central site r = 1 and moving to the left 
{Vp < 0) pass the site r two times (once before the reflection at the left boundary and 
once after the reflection) before they return (now with Vp > 0) to their starting point. 
These two times are ti = \r\/vp and t2 = Tp — ti. Therefore the probability that this 
QP trajectory passes r at times t < Tp exactly once is 

t<ti 

1 ti<t<t2 (34) 

t2<t<Tp 

The associated partner (of the QP pair) that moves to the right {vp > 0) passes the site 
r only after reflection at the right boundary and returns to the starting point, at times 
= Tp + ti and t^ = 2Tp — ti. Therefore the probability that this QP trajectory passes 
r at times t < 2Tp exactly once is 

' t<t3 

1 ts < t < t4 (35) 
ti<t<2Tp 

For t > 2Tp one makes use of the Tp-periodicity of qp{t): 

qp{t + 2nTp) = qp{t), (n = l,2,...). (36) 

With ( l32l) . ( l33l) . (1341I36P one obtains mi{t) via numerical integration (or summation over 
the discrete values for a lattice of finite size L), where we take the exact profile before 
the quench (m;(t = 0)). 

In Fig. [Ulthe comparison of this QP calculation with our exact data is shown. Fig. 
[TTK displays the dynamical evolution of the local magnetization ai h = 0.2 close to the 
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Figure 11. (Color online) Local order parameter mi{t) within the FM phase {h < 1) 
and comparison with QP calculation. A: mr{t) at h — 0.2 as a function of the distance 
r from the defect for different fixed times t. B: mi{t) at h = 0.5 as a function of I for 
different fixed times (L = 256, the defect is at Z — 128). The full lines are the exact 
data and the points are the predictions of the QP calculation. C-D: mi{t) for different 
fields /i < 1 as a function of the rescaled time t/T{h), where T{h) — L/v^iix{h) — L/h, 
i.e. T{h = 0.5) = 512, T{h = 0.2) = 2560, and T{h = 0.1) = 5120. The exact data and 
the QP prediction can be discriminated by the smoother behavior of the latter. 



defect. The system size, L = 256, is sufficiently large such that the boundary effects are 
not yet visible for the times shown and the data are representative for an infinite system. 
One sees that the initial magnetization drop at the defect (the defect spins are surface 
spins at t = 0) is quickly ffiled and the profile approaches a constant magnetization 
nib = ^i^L,t^oofnr{t)- Since the probability qp{t) for a QP with momentum p to pass 
any site in the bulk is 1 in the the limit L — > oo and t — > oo the bulk magnetization is 
predicted by the QP picture, according to f l32l) and (l3T]) . to be 

mfo = mfeexp(-l/0 , (37) 

with 

1 _ 1 

^ is identical with the nonequilibrium correlation length measured during global quench, 
as defined in Eq.([n]). For fp{h) as in Eq. fl33p . which is exact for small h it is given by: 
1/^ = h'^/A. For larger values of h corrections to the kink-like excitations should by 
summed, which leads to the value ^ [2D]: 1/^ = - In [(l + Vl - h^^ /2] . For h = 0.2 



dpfpih) . (38) 
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one gets = 0.985 which agrees well with the approximately constant magnetization 
at t = 640 in Fig. [TTK. We checked that also for larger values of h the QP prediction 
for the asymptotic bulk magnetization nib in eq.( l37|) is good. 

Fig. [TTB compares the exact data for the local magnetization profile of a finite 
system for different times after the quench with those of the QP prediction at h = 0.5. 
The agreement is very good for the times shown and a few sites away from the defect 
(r = |L/2 — /| > 5). The rapid filling of the initial sharp dip of the profile at the center 
(see Fig. [TT1A-) is not correctly captured by the present QP picture. To understand this 
process one should assume, that the QP-s are emitted in a region of size C,s around 
the defect and these QP-s are quantum entangled and these correlated particles are 
responsible for the rapid increase of the magnetization at the defect. 

The predicted (quasi) -period is governed by the maximum QP velocity and given 

by 

T{h) = L/v^^, ^ L/h . (39) 

This agrees well with the dynamical behavior of m/(t) shown in[TT]C-D displaying the 
exact data and the QP prediction for m/(t) at fixed sites / for different fields h < 1 
as a function of the rescaled time t/T{h). The agreement for the sites shown is good, 
at longer times close to the defect (/ ^ L) deviations occur (see Fig. [TTD) due to the 
mechanism described above. 



5. Quench in the paramagnetic phase 

In the paramagnetic phase {h > 1) the magnetization is vanishing in the thermodynamic 
limit, however, in a finite system there is a finite, L-dependent magnetization. Its 
temporal evolution for /i = 2 is depicted in Fig. [12] for L = 128. In contrast to the 
dynamics at the critical point (c.f. Fig. HI) the evolution of mi{t) is not periodic in the 
paramagnetic phase and approaches a stationary profile characteristic for the system 
without defect. 

Our results for the spatial correlation function Cj(ri,r2) within the paramagnetic 
phase are depicted in Fig. [131 We find that for fixed ri and r2 an asymptotic power law 
in t with an exponent close to 3/2: 

Ctiri,r2) oc t"^/^ for t > ri, jra] . (40) 

Finally the autocorrelation function after a local quench in the paramagnetic phase 
decays as in the ground state without a defect, i.e. 

Gt{r) oc t-'/' . (41) 

see Fig. [H 



Non- equilibrium quantum dynamics after local quenches 



19 




Figure 12. (Color online) Temporal evolution of the local magnetization profiles 
after a local quench within the disordered (PM) phase, here h = 2, L = 128. 
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Figure 13. (Color online) Spatial correlation function Ct(ri,r2) as a function of t 
for local quenches in the paramagnetic phase (ft, > 1) in a log-log plot. The straight 
lines correspond to the power law A: Different pairs of (ri,r2) with the same 

spatial difference at ft = 3 (L is 256). B: ft = 2, i = 256. C: ft = 5, i = 256. D: 
ft = 5, L = 460. 
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Figure 14. (Color online) Autocorrelation function Gr{t) after a quench in the 
paramagnetic phase for r = 10 (i.e. for a site in a distance 10 from the defect) for 
different values of h in comparison with Gr{t) for a system without defect ("nd"). 



6. Conclusions 

We have studied tlie temporal evolution of different observables in the transverse Ising 
chain following a local quench: for t < the system consisted of two disconnected 
halves which are joined together for t > with the uniform bulk coupling, 7^/2 = 1- 
We have measured the magnetization profile, mr{t), as well as the correlation, Ct{ri, r2) 
and the autocorrelation function, Gr{ti,t2). We have concentrated on the properties 
of local quench in the critical state, but some calculations are also performed in the 
ferromagnetic and in the paramagnetic phases, too. 

For critical local quench several conjectures about mr{t) and Ct{ri,r2) are known 
through conformal field theory, which are valid in the thermodynamic limit and in the 
continuum approximation. Our exact finite lattice results have confirmed the conformal 
conjectures, except for the early time behavior of the correlation function in which 
the two reference points are at different sides of the defect. We have also studied 
systematically the finite-size effects, in particular we have made conjectures about the 
form of the magnetization profiles, both in space Eq.f llSp and time Eq. (fT7|) . which 
probably can be derived in some way, e.g. through conformal field theory. 

Our results are explained within the frame of a quasi-particle picture, in which 
during the quench kink-like excitations are created at the defect, which move semi- 
classically, with a momentum-dependent velocity and result in the reduction of the 
order-parameter in the system. In the ferromagnetic phase this type of semi-classical 
calculation has lead to exact results, at least for a small transverse field, h. We expect. 
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however, that following the same method as in the case of global quench, one can sum 
the higher order contributions and in this way one obtains exact results about the 
stationary value of the magnetization profile for general value of h < 1. This involves 
the non-equilibrium correlation length, ^, as measured at a global quench, see Eq. (l37|) . 

Our investigations can be extended and generalized in different directions. First, 
one can consider another quantum spin chains, for which the conformal conjectures at 
the critical point are expected to be satisfied in the same way as for the transverse Ising 
chain. In the ferromagnetic phase the relation in Eq.( l37|) is expected to be valid and 
in this way one can measure the non-equilibrium correlation length in an independent 
procedure. A second way to generalize our results is to use different forms of the local 
quench. One possibility is to use a non-zero defect coupling between the two subsystems 
in the initial state and/or to have a non-uniform defect coupling, 7^ = 1, in the 
final state. In the transverse Ising chain the local critical exponents are continuous 
function of the strength of the defect |^ HU |l5], so that Xs for a decoupled system 
should be replaced by a defect exponent, Xd- Finally, we can also study local defects 
with a more complicated structure, which involve several lattice sites. 
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